A HEURISTIC ALGORITHM FOR COMPUTING THE 
POINCARE SERIES OF THE INVARIANTS OF BINARY FORMS 



DRAGOMIR Z. BOKOVIC 

Abstract. We propose a heuristic algorithm for fast computation of the 
Poincare series P n (t) of the invariants of binary forms of degree n, viewed 
as rational functions. The algorithm is based on certain polynomial identities 
which remain to be proved rigorously. By using it, we have computed the 
Pn(t) for n < 30. Dixmier proposed long ago three conjectures concerning 
these Poincare series. We verify that the first two of them are valid in the 
above range. As a supplement to this note, we provide a table from which one 
can easily construct the functions P n (t) for 3 < n < 30. Only the portion of 
the table covering the range 3 < n < 20 is actually displayed in the paper. 



1. Introduction 

Let V(V) be the algebra of polynomial functions on a finite-dimensional complex 
vector space V. Let V n (V) C V be the space of homogeneous polynomials of degree 
n. We shall assume from now on that V = C and refer to V n (C 2 ) as the space of 
binary forms of degree n. By fixing a basis, we shall view V(C 2 ) as the polynomial 
algebra C[x,y] in the two coordinate functions x and y. The natural action of 
GL2(C) on V induces an action on V(C 2 ) such that each subspace V n (C 2 ) is 
stable and the representation of GL2(C) on it is irreducible. 

To simplify the notation, we set G — SLi2(C). One of the main themes of the 
classical invariant theory was the study of the algebra of G- invariants of V(V n {C 2 )), 
n = 1,2,3,... Denote this subalgebra by T„ = V{V n {C 2 )) G . This notation is 
borrowed from the article of R. Howe [5], which gives a modern point of view on 
this classical subject. The Poincare series (also known as the Hilbert series) of I n 
is the formal power series 

P„(t)=]Tdim(^)t fc , 

fc>0 

where 

i k n = i n nv k (P n (c 2 )). 

It is well known that this series is in fact the Taylor series of a rational func- 
tion, which we denote also by P n {t). We consider here the problem of computing 
efficiently these rational functions. 

There are two methods that we find in the recent literature on this subject. First 
there is a method described in Howe's article and implemented on a computer by 
P. Saly, Jr. It is based on the classical formulae, due to Cayley and Sylvester, for 
the coefficients of the above Taylor series and the fact that the concrete form for 
the denominator of P n (t) is known. Indeed, Dixmier's Conjecture 1 in gives 
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simple expressions for the lowest term denominators of the P n 's. (See also Section 
El below.) 

The second method is based on an explicit but complicated formula for P n (t) 
discovered by T. Springer This formula has been used by Brouwer and Cohen 
PP and also by Littelmann and Procesi [Hj. 

Our method is completely different and is based on certain conjectural poly- 
nomial identities (which are not explicitly given). In spite of this vagueness, the 
conjectural formulae can be used to write an efficient algorithm to compute the 
P n 's for small n's. The algorithm saves time by avoiding the usual tedious proce- 
dure of computing many residues one at the time. We used Maple [7] to obtain the 
formulae for P n for n < 30. The memory requirements are very modest and the 
computations can be performed on a PC. 

We did not have access to the tables of Brouwer and Cohen (for n < 17) to 
compare them with our tables. At the end of his article Howe gives two concrete 
coefficients of t k in the numerator of P n (t), one for n = 19 and k — 160 and the 
other for n = 20 and k — 84. The denominators are specified as n 7 =2(^ — ^"O ^ or 



n = 19, and ^=2(1 ~~ ^) f° r n = ^0. We find that the first coefficient is correct 
but the second one should be multiplied by 2. 

In a separate file, we provide the list of the P„'s for n < 30 as they may be of 
interest to other researchers. (I wish I had access to such a list when I was doing 
these computations.) More details about this list are given in Sectional 

2. The main actors 

Let z and t be independent commuting variables and Z[z,<] the corresponding 
polynomial ring with integer coefficients. 

Let n > 3 be an integer. Set n = 2s — 1 if n is odd and n = 2s if n is even. The 
main actors in this paper are the three polynomials p n ,q n £ Z[z,t] and r n 6 Z[i]. 
The first two are defined, for odd n, by 

s s 

p n (z, t) = H(l - tz 2 ^ 1 ), q n (z, t) = Y[(z 2i l - t), 
and, for even n, by 

s s 

Pn(z,t) = - tz 2i ), q n (z,t) - l[(z 2i - t), 

i=l i=l 

In both cases, s is the t-degree of p n (and q n ), i.e., its degree as a polynomial in 
the variable t. Denote by m the z-degree of p n (and q n ). Thus, m = s 2 if n is odd 
and m = s(s + 1) if n is even. 

Observe that 

(2.1) z m p n {z-\t) = q n (z,t), {-typ^z^-^^q^z^t), 

and, consequently, 

z m (-t) s p n (z-\r 1 )= Pn (z,t). 

These formulae remain valid when the letters p and q are interchanged. 
The third polynomial is defined by 

n-l 

r n (t) = l[(l-t 2i ) 
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if n is odd, and by 

n-l 

r n (t) = (l + i)n(l-0 

if n is even. 
Finally, set 

Vn {z,t) = z m - 2 (z 2 -l)r„(t). 

It is easy to verify that 

(2.2) 2 2m -Vn(*-\*) = -¥>n(z,*), ^ n (z,t- 1 ) = (-l)Vn(«,t). 

and, consequently, 

z^-VV^z- 1 ,;- 1 ) = (-1)*+Vn(«,t), 

where d is the degree of r n (t). Thus d = 2s(n — 2) for n odd and d — s(n — 1) for 
n even. 

3. Basic conjecture 

In this section we state our basic working conjecture, which will be taken for 
granted in the remaining part of the paper. It is supported by extensive computer 
calculations. 

Let 7„ =< p n , q n > be the ideal of Z[z, t] generated by p n and q n . We conjecture 
that <p n G I n . More precisely, 

Conjecture 3.1. There exist unique polyjiomidls b n G Z[z, t] of z-degree m — 2 

(3.1) <Pn{z, t) = a„(z, t)q n (z, t) + b n (z, t)p n (z, t). 

(The uniqueness is clear because p n and q n are relatively prime.) 

Let us give an example. If n = 5 then s = 3, to = 9 and (p = aq + bp, where 

p = zV-i)(i-; 4 )(i-t 6 )(i-; 8 ), 

p = (l-iz)(l-tz 3 )(l-iz 5 ), 
q = (z-t)(z 3 -t)(z 5 -t), 

a = {l-t 6 + t 12 )-t(l + t 2 -t 6 )(l-t 2 -t 6 )z + t 2 (l-t 6 +t 8 )z 2 

-t(l - i 4 )(l - t 6 - i 8 )z 3 + t 2 (l - t 4 )(l + t 2 - i 8 )z 4 - t 5 (l -t 2 + t 8 )z 5 

+t 2 (l + t 4 - t 6 )(l -t 4 - t 6 )z 6 - t 3 (l -t 6 + t 12 )z 7 , 
b = t 3 (l-t 6 + t 12 )-t 2 (l + t A -t 6 )(l-t 4 -t 6 )z + t 5 (l-t 2 +t 8 )z 2 

-t 2 (l - t 4 )(l +t 2 - t 8 )z 3 + t(l - t 4 )(l -t 6 - t s )z A - t 2 (l -t 6 + t s )z 5 

+t(l + t 2 - t 6 ){l -t 2 - t 6 )z 6 - (1 - t 6 + t 12 )z 7 . 

By applying the substitution z — ► z -1 to l|3.1|l and by using the identities (|2.1|) 
and H2.2f) . we obtain that 

(3.2) b n {z,t) = -z m - 2 a n {z-\t). 
Similarly, the substitution t — > t^ 1 gives 

(3.3) b n {z,t) = (-l) n+ '1^'a n (z,t- 1 ). 
From the last two equations we deduce that 

(3.4) z m ~ H d - s a n {z-\ r 1 ) = (-l)" +s+1 a„(z, t). 
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We can write 

m— 2 m— 2 

(3.5) a n (z,t) = a k {t)z k , b n {z,t) = £ (3 k (t)z k , 

2 = i = 

where ctk,f3k G Z[t]. Then the equation 1)3. 2[1 gives 

(3.6) /8*(*) = -Om-a-fcW, 0<fc<m-2 
and 1|3.4[1 gives 

(3.7) a„ l _ 2 _ fc (i) = (-l)"+ s+1 < !i - s a fe (<- 1 ), < k < m - 2. 
Hence, the equation l|3.1[l can now be rewritten as 

m-2 

(3.8) <p n (z,t) = a k(t) {z k 1n{z,t) - z m - 2 - k Pn (z,t)) . 

k=0 

By setting z = in (|3.8(l . we obtain that a m -2(£) = (— £) s ao(£). By combining 
this with the equation 1)3. 7)1 for fc = 0, we obtain that 

(3.9) ao(t) - [-l) n+1 t d - 2 'ao(t- 1 ). 

Apparently the computation of the polynomials a n (z,t) is a difficult job. For- 
tunately, we do not need to know these polynomials explicitly but only the poly- 
nomials a„(0, t) = ao(t). The main point of our algorithm is that the polynomials 
ceo(t) can be computed efficiently. 

For the curious reader, let us throw in one more intriguing conjecture (which 
plays no role whatsoever in this paper). 

Conjecture 3.2. J n flZ[t] is the principal ideal ofZ[t] generated by the polynomial 

n—1 n — 1 n—1 

i—1 i—1 i—1 

according to whether n is odd, congruent to 2 modulo 4, or divisible by 4. 
This conjecture makes sense (and is true) also when n is equal to 1 or 2. 

4. The integral formula 

The Molien-Weyl integral formula for P n (t) can be written in the following form 
(see EJ) 



2™ J M=1 z 

fn(z,t) 



'1*1=1 

where 

nLo(i-te"- 2fe )' 

The integration is to be performed over the unit circle in the counterclockwise 
direction and it is assumed that \t\ < 1. 
Let us define yet another function 

z m- 2(z 2 _ ^ 



»(*,*) = 



Pn(z,t)q n (z,t)' 
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For odd n, f n (z,t) — g n (z,t) and, for even n, f n (z,t) = g n (z,t){\ ~t)~ l . Since the 
integration variable is z, the factor (1 — i) -1 can be inserted at the very end of the 
computation. Thus it suffices to compute the integral 

_ , . 1 f dz 

Qn(t) = 7T~ / 9n(Z,t) . 

27r * J\z\ = l z 



The formula gives 

r n (t)g n (z,t) = 

Hence, 

r n (t)Q n (t) = 



i(z,t) b n (z,t) 



p n (z,t) q n (z,t) 
1 f a n (z,t) dz 1 f b n {z,t)dz 



2tti J\ z \ = i Pn(z,t) z 2ni J\ z \ =1 q n {z,t) Z 

The second integral is since all the poles are inside the unit circle (and the residue 
at co is 0). The integrand of the first integral has all of its poles outside the unit 
circle except for the simple pole at z = 0. Since a n (0, i) = ao(t) and p n (0,t) = 1, 
we obtain that 

(4.1) 

Hence, the computation of P n (t) is reduced to computing the polynomial ag (t) . 
Since Q„(0) = r„(0) = 1, we deduce that ao(0) = 1. The formula l|3.9|l now implies 
that ao has degree d — 2s. This is in agreement with the well known fact (it was 
known to Sylvester [Hj) that the rational function P n (t) has degree — n — 1. 

If n is odd, then ao(t) is a palindromic polynomial, i.e., its coefficients are sym- 
metric about the midpoint. If n is even, then ao(l) — 0, i.e., ao(t) is divisible by 
1-t. 

Let us emphasize that the main feature of our algorithm described above is that 
we compute the above integral at one fell swoop, avoiding the tedious procedure of 
computing all the residues for the poles inside the unit circle one at the time. 

5. HOW TO COMPUTE a Q ? 

Let me start by quoting the Rule #1 from the list of useful tips on Jean-Charles 
Faugere's home page: "In a first time try to compute modulo p where p is a small 
prime." I learned this rule on my own (i.e., the hard way). 

To compute ao(t) we used the Maple package called "Linear Algebra" and its 
subpackage called "Modular". First we choose a big prime, I, in the data type 
integer [4]. The largest one is I — 65521. We perform computations modulo this 
prime. We generate random integer mod I inputs for the variable z and plug it into 
our equation l|3.8l) . We need only m — 1 inputs since there are m — 1 unknown a^'s. 

Thus we obtain a system of m — 1 equations which arc linear in the a^'s. But 
these are still polynomial equations as they contain the variable t. We chose a 
random value, say t modulo Z, for the variable t and plug it into this system of 
equations. Now we have just a small system of linear equations to solve for the 
unknowns afc(r). This job can be easily handled by "Modular". We are only 
interested in the value cco(t). 

Next recall that we know the degree of ao; it is equal to d — 2s. Thus we have to 
repeat the above calculation d — 2s + 1 times to get the values ao{r) for d — 2s + 1 
different inputs r. Having done this, the d— 2s + 1 unknown coefficients of ao(t) can 
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be easily computed by solving the corresponding Van der Monde system of linear 
equations. Of course, this means that at this stage we know these coefficients 
modulo our prime I. 

Next we have to repeat the whole calculation above with several distinct primes. 
We used the largest seven primes available: 

6552f , 655f9, 65497, 65479, 65449, 65447, 65437. 

This is necessary when n is near the high end of our range because some of the 
coefficients of oio(t) for n = 29 have 30 digits. The remaining task is to lift these 
modular solutions to a genuine solution over Z. For this purpose, "Modular" pro- 
vides the "ChineseRemainder" command which makes the task very easy. 

In the hardest case, n = 29, the computation of ao(t) modulo the above seven 
primes (performed one after the other) took almost 5 hours on the SunBlade work- 
station running a single R10000 CPU at 250 MHz with 8 GB of RAM. 

Certainly, if Coniecture l3. II was false these computations, which involve so many 
random choices, would not produce a meaningful result. Moreover, the polynomials 
ao(t) have certain symmetry properties. As mentioned earlier, if n is odd then ao(t) 
is a palindromic polynomial. If n is even, then ao(^) is divisible by I — t, according 
to one of the Dixmier's conjectures. Finally, if n is divisible by 4, then ao(t) is also 
divisible by I + 1 3-1 , again by Dixmier. This is all true in the cases that we have 
computed. 

If n > 3 then t = I is a pole of P n (t) of order n — 2 and the limit of (1 — t) n ~ 2 P n (t) 
as t — > I has been computed by Hilbert 0]. We have verified that our formulae for 
P n {t) agree with his result. 

6. Description of the coefficient table 

There exist unique relatively prime polynomials A n ,B n G Z[t], with A n (Q) = 
B n (0) = 1, such that 

In this notation, Dixmier's Conjecture I asserts that B n (t) — r n (t)/(l+t s ) if n is 
divisible by 4, and B n (t) = r n (t) otherwise. He verified this conjecture for n < 16 
and our computations show that it is valid for n < 30. 
We propose the following additional conjecture: 

Conjecture 6.1. The polynomials A n (t) are irreducible for n > 5. 

Since the denominators B n (t) are easy to write down, we provide only the infor- 
mation necessary to construct the numerators A n (t) for 3 < n < 30. 

If n is odd, then A n (t) = oto(t). Now assume that n is even. Then ao(t) = 
(1 - i)ao(f) for some a G Z[t] and P n (t) = a Q (t)/r n (t). Thus A n (t) = a (t) if 
n = 2 (mod 4). Finally, if n = (mod 4), then ao(t) is divisible by 1 + and 
we have A n (t) = a (t)/(l + t 3 - 1 ). 

The degree of A n is always even and we shall denote it by 28. We have 5 = 
2s(s - 2) for n odd, 25 = s(2s - 3) - 1 for n = 2 (mod 4), and 5 = s(s - 2) for 
n = (mod 4). These numerators are palindromic polynomials 

A n {t) = c + cit + ■ ■ ■ + cs-it*- 1 + c s t s + c S -it s+1 + ■■■ + erf 25 ' 1 + c a t 28 

and so it suffices to store roughly half of their coefficients. 
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For each n in the range 3 < n < 30 we record the coefficients cq — 1, c\, . . . , eg 
of A n in that order. They are stored as a Maple table A[3..30]. 

For example, if n — 5 then s — 3, 6 — 6, A$ = 1 — t 6 + t 12 . This is recorded in 
our table as 

A[fi] - [1,0,0,0,0,0,-1]. 
Assume that n > 3 is odd. Then it is known that ci — for odd «'s as well as for 
i = 2. Dixmier 3 shows that C4 and eg are nonnegative, and so are cq for n > 15 
and C10 for n > 9. We can prove that in fact 

n — 3 



It may be useful to obtain similar formulae for ce, c%, etc. 

7. Appendix: Table A[3..20] 

We display here the coefficients cq = 1, ci, . . . , ca of the numerator A n (t). Recall 
that the degree of A n (t) is 26. The display includes only the values n = 3, 4, . . . , 20. 
The full table with the range n = 3, 4, . . . , 30 is available from the author on request. 
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A[3] = [1] 
A[4] = [1] 

A[5] = [1,0,0,0,0,0,-1] 
A[6] = [1,1,0,-1,-1] 

A[7] = [1, 0, 0, 0, 0, 0, -1, 0, 2, 0, -1, 0, 5, 0, 2, 0, 6] 
A[8] = [1,1,0,-1,-1,0,1,1,1] 

A[9] = [1, 0, 0, 0, 1, 0, -1, 0, 5, 0, 3, 0, 18, 0, 15, 0, 44, 0, 43, 0, 82, 0, 76, 0, 122, 0, 
107,0,147,0,119] 

A[10] = [1, 1, 0, -1, -1, -1, 2, 2, 4, 6, 10, 10, 16, 15, 16, 16, 19, 15] 
A[U] = [1, 0, 0, 0, 1, 0, -1, 0, 10, 0, 11, 0, 58, 0, 85, 0, 222, 0, 336, 0, 660, 0, 951, 0, 
1589, 0, 2154, 0, 3188, 0, 4080, 0, 5510, 0, 6633, 0, 8310, 0, 9443, 0, 11059, 0, 
11894, 0, 13094, 0, 13319, 0, 13852] 

A[12] = [1, 1, 0, 0, 1, 1, 3, 7, 11, 15, 25, 35, 50, 70, 95, 121, 156, 191, 229, 266, 304, 
332,358,371,377] 

A[13] = [1, 0, 0, 0, 1, 0, -1, 0, 19, 0, 31, 0, 157, 0, 321, 0, 885, 0, 1756, 0, 3794, 0, 
6856, 0, 12788, 0, 21324, 0, 35633, 0, 55326, 0, 85174, 0, 124064, 0, 178645, 0, 
246238, 0, 334814, 0, 439321, 0, 568305, 0, 712862, 0, 881834, 0, 1061455, 0, 
1259989, 0, 1459221, 0, 1666984, 0, 1860904, 0, 2049854, 0, 2209072, 0, 2349306, 
0, 2446352, 0, 2514111, 0, 2530530] 

A[U] = [1, 1, 0, -1, 0, 0, 5, 8, 21, 34, 69, 107, 193, 295, 477, 705, 1064, 1490, 2121, 
2865, 3876, 5040, 6535, 8186, 10217, 12379, 14898, 17490, 20381, 23197, 26216, 
28992, 31799, 34192, 36461, 38127, 39559, 40263, 40636] 
A[U] = [1, 0, 0, 0, 2, 0, 0, 0, 32, 0, 76, 0, 378, 0, 995, 0, 3048, 0, 7294, 0, 17681, 0, 
37736, 0, 78903, 0, 152321, 0, 285968, 0, 507762, 0, 876759, 0, 1451423, 0, 2341739, 
0, 3653241, 0, 5568497, 0, 8254649, 0, 11983447, 0, 16987847, 0, 23631274, 0, 
32196429, 0, 43116834, 0, 56681420, 0, 73342055, 0, 93320393, 0, 117007543, 0, 
144461993, 0, 175919353, 0, 211175615, 0, 250222591, 0, 292516508, 0, 
337751801, 0, 385016863, 0, 433713649, 0, 482605505, 0, 530877973, 0, 
577086324, 0, 620343376, 0, 659172312, 0, 692798202, 0, 719914717, 0, 
740045690, 0, 752239053, 0, 756462172] 

A[16] = [1, 1, 0, 0, 1, 3, 10, 18, 38, 79, 153, 278, 514, 891, 1523, 2528, 4072, 6367, 
9772, 14572, 21306, 30498, 42785, 58863, 79666, 105852, 138459, 178275, 226114, 
282612, 348514, 423785, 508764, 603037, 705993, 816528, 933532, 1054746, 
1178503, 1302160, 1423129, 1538529, 1645908, 1741955, 1824693, 1891546, 
1940709, 1970683, 1980976] 
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A[17] = [1,0, 0, 0, 2, 0, 0, 0, 50, 0, 156, 0, 844, 0, 2716, 0, 9280, 0, 26055, 0, 
70846, 0, 173224, 0, 405183, 0, 883551, 0, 1847356, 0, 3669433, 0, 7024773, 0, 
12919848, 0, 23019526, 0, 39697193, 0, 66608244, 0, 108748704, 0, 173371011, 
0, 270001994, 0, 411791616, 0, 615371715, 0, 902700319, 0, 1300556398, 0, 
1842885348, 0, 2569619659, 0, 3529481145, 0, 4777707107, 0, 6379225544, 0, 
8404807944, 0, 10934524315, 0, 14051849433, 0, 17847385164, 0, 22410522390, 
0, 27833575972, 0, 34200807232, 0, 41593316021, 0, 50075498973, 0, 
59701312756, 0, 70498510402, 0, 82477120228, 0, 95612069584, 0, 
109854817368, 0, 125115124389, 0, 141276752124, 0, 158179073390, 0, 
175637520454, 0, 193425539961, 0, 211299997259, 0, 228983624510, 0, 
246195050087, 0, 262631400199, 0, 278002047047, 0, 292010632283, 0, 
304391087240, 0, 314888097678, 0, 323293686822, 0, 329425050810, 0, 
333160013567, 0, 334411422423] 

A[18] = [1, 1, 0, -1, 1, 2, 11, 20, 60, 122, 292, 573, 1199, 2264, 4307, 7692, 13639, 
23121, 38688, 62619, 99808, 154969, 236962, 354532, 522744, 756614, 1080091, 
1517149, 2103843, 2875718, 3883971, 5178127, 6826649, 8894132, 11466985, 
14623600, 18466440, 23083747, 28588450, 35071007, 42645944, 51394256, 
61419462, 72778864, 85549508, 99748823, 115411657, 132499497, 150990759, 
170778600, 191774520, 213796608, 236686776, 260188811, 284080692, 
308043009, 331804545, 355002792, 377341521, 398442572, 418014893, 
435697784, 451237581, 464324429, 474773647, 482354991, 486977592, 
488512665] 
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A[19] = [1, 0, 0, 0, 2, 0, 1, 0, 76, 0, 296, 0, 1763, 0, 6738, 0, 25712, 0, 82893, 0, 
252012, 0, 694765, 0, 1807368, 0, 4392969, 0, 10154779, 0, 22296093, 0, 
46930799, 0, 94802787, 0, 184822672, 0, 348299741, 0, 636837951, 0, 1131559371, 
0, 1959027689, 0, 3309329549, 0, 5465457626, 0, 8835569146, 0, 14002953513, 0, 
21778902561, 0, 33281420196, 0, 50015102190, 0, 73986433683, 0, 107815040622, 
0, 154890725357, 0, 219515780752, 0, 307103168571, 0, 424341546490, 0, 
579425240340, 0, 782222455286, 0, 1044516156680, 0, 1380143077777, 0, 
1805211986577, 0, 2338168202449, 0, 2999945060393, 0, 3813898895206, 0, 
4805843429850, 0, 6003797685274, 0, 7437855144286, 0, 9139694108669, 0, 
11142247265979, 0, 13478942963820, 0, 16183166152183, 0, 19287230767266, 
0, 22821658964018, 0, 26813925195712, 0, 31287624164380, 0, 36261071290807, 
0, 41746472475910, 0, 47748507080547, 0, 54263646677747, 0, 61278874674037, 
0, 68771319597259, 0, 76707272862777, 0, 85042295486194, 0, 93720678352722, 
0, 102676146969267, 0, 111831862933442, 0, 121101792448393, 0, 
130391278393710, 0, 139599058259333, 0, 148618354053166, 0, 
157339436459345, 0, 165651094224795, 0, 173443559767368, 0, 
180610146679903, 0, 187050288671933, 0, 192671089114197, 0, 
197390193869982, 0, 201136986494180, 0, 203855025614692, 
0, 205502649265697, 0, 206054755643582] 

A[20] = [1, 1, 0, 0, 2, 5, 17, 40, 100, 232, 544, 1199, 2599, 5365, 10770, 20867, 
39312, 71826, 128004, 222286, 377375, 626606, 1019690, 1627231, 2550571, 
3929105, 5955664, 8888409, 13072800, 18958777, 27131530, 38333433, 53503793, 
73805811, 100673389, 135839979, 181391167, 239789025, 313924970, 407134237, 
523239787, 666546890, 841871660, 1054502755, 1310204347, 1615137665, 
1975827494, 2399034646, 2891684840, 3460685330, 4112816074, 4854497808, 
5691649536, 6629418196, 7672023583, 8822464074, 10082371655, 11451725459, 
12928740940, 14509620528, 16188514360, 17957332700, 19805806095, 
21721392176, 23689449970, 25693256426, 27714311797, 29732461463, 
31726325603, 33673521588, 35551188212, 37336279422, 39006154571, 
40538894181, 41913908026, 43112232559, 44117095575, 44914144829, 
45491928155, 45842000355, 45959277535] 
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